Floquet spin states in graphene under ac driven spin-orbit interaction 
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We study the role of periodically driven time-dependent Rashba spin-orbit coupling (RSOC) on 
a monolayer graphene sample. After recasting the originally 4x4 system of dynamical equations as 
two time-reversal related two-level problems, the quasi-energy spectrum and the related dynamics 
are investigated via various techniques and approximations. In the static case the system is a 
gapped at the Dirac point. The rotating wave approximation (RWA) applied to the driven system 
unphysically preserves this feature, while the Magnus-Floquet approach as well as a numerically 
exact evaluation of the Floquet equation show that this gap is dynamically closed. In addition, 
a sizable oscillating pattern of the out-of-plane spin polarization is found in the driven case for 
states which completely unpolarized in the static limit. Evaluation of the autocorrelation function 
shows that the original uniform interference pattern corresponding to time-independent RSOC gets 
distorted. The resulting structure can be qualitatively explained as a consequence of the transitions 
induced by the ac driving among the static eigenstates, i.e., these transitions modulate the relative 
phases that add up to give the quantum revivals of the autocorrelation function. Contrary to the 
static case, in the driven scenario, quantum revivals (suppresions) are correlated to spin up (down) 
phases. 

PACS numbers: 81.05.ue, 71.70.Ej, 72.25.Pn 



I. INTRODUCTION 

One of the key features of relativistic (massless) free 
particle states is that they evolve, at least in effectively 
one-dimensional situations, in time without spreading. 
This in turn relies on the linear nature of the relativis- 
tic dispersion relation which for photons reads u>(k) = 
ck, with c the speed of light. The condensed mat- 
ter relativistic-particle analog is found in the low en- 
ergy approximation (long wavelength) of single layer 
graphene where the chiral massless particles move with 
a speed Vp ss c/300 1 . This linear spectrum provides 
graphene with remarkable transport properties such as 
high mob ility ^, Klein tunneling^ and unconventional spin 
Hall effecpEl. The later stems from the interplay between 
intrinsic spin-orbit interaction and the coupling extrinsi- 
cally induced by external gate voltages or an appropriate 
substrate. This extrinsic couplingSHSl, so called Rashba 
spin-orbit (RSOC) has been als o fou nd to give rise to 
spin polarizatiorP and relaxation ! 10 * 11 ! effects. 

Although the role of static RSOC on graphene has been 
extensively discussed in the literature, to the best of our 
knowledge, the role of periodically driven time dependent 
RSOC on graphene samples has not been analysed so far. 
Yet, recent works have focused on the dynamical features 
of charge currents induced by means of time dependent 
extrinsic spin-orbit interaction on mesoscopic semicon- 
ductor quantum rings where Rabi oscillations are shown 
to appear as well as collapse and revival phenomena^. 
The main motivation of our work is twofold: First at 
all, we are interested in determining the feasibility of ac 
driven fields to generate and modulate a finite spin polar- 
ization of carriers in graphene for states that under static 
conditions remain unpolarized. In addition, we are also 



interested in the dynamical modulation of the effective 
Rashba coupling strength A = \/Ml which would allow 
to explore regimes beyond the static limit domain. 

Taking advantage of the periodicity of the problem, 
the evolution equations can be solved via Floquet the- 
ory. A standard approach here consists of expressing the 
Hamiltonian in a Fourier mode expansion leading to an 
infinite-dime nsiona l eigenvalue problem for the so-called 
quasi-energie^^. This quasi-energy spectrum carries 
nontrivial information on the topological nature of the 
system under study*^, and for semiconductor quantum 
wells with a zincblende structure it has been recently 
shown that ac driving can induce a topological phase 
transition^. 

Practically, in order to treat the infinite eigenvalue 
problem one has to truncate at an order of the harmonic 
expansion chosen appropriately to yield well-converged 
results. An alternative approach to the Floquet prob- 
lem which does not rely on Fourier expansions has been 
devised by Magnus^. This method appears to be some- 
what less popular and amounts in formulating the time 
evolution operator as the exponential of a series of nested 
commutators. It has the virtue of both preserving uni- 
tarity at any order in the series expansion (in contrast 
to truncation of the Dyson series within a perturbative 
approach) and avoids the infinite-dimensional eigenvalue 
problems. Following^ we will make use of the Magnus 
expansion approach combined with Floquet theory in or- 
der to generate semi-analytical solutions of the dynamics 
induced by periodic RSOC. 

Since RSOC couples the spin and pseudospin degrees 
of freedom the problem is, for a given wave vector, gener- 
ically four-dimensional. However, by an appropriate uni- 
tary transformation, the evolution equations can be re- 
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cast as a set of two equivalent two-level Schrodinger equa- 
tions related by time reversal. In this way we can explic- 
itly analyse which static states get dynamically coupled. 

Our main results are the following: The ac driven 
RSOC induces a quasi-energy spectrum where the origi- 
nal gap due to static spin-orbit coupling is dynamically 
closed. In particular, at the Dirac point (k = 0) the 
dynamics is exactly solvable with zero quasi-energy Flo- 
quet states. This quasi-energy spectrum is two fold- 
degenerate as a consequence of time reversal invariance of 
spin-orbit interaction, and the closing of the original gap 
is due to the destructive interference induced among the 
initially uncoupled positive and negative energy RSOC 
eigenstates. Then we show that sizeable alternating out 
of plane spin polarization ensues on states that under 
static conditions remain unpolarized. We also find that 
the uniform interference pattern shown by the autocorre- 
lation function for static RSOC gets distorted due to the 
interlevel mixing of the static eigenstates which dynam- 
ically modulates the relative phases that add up in the 
quantum revivals of the autocorrelation function. In the 
driven case quantum revivals (suppressions) are directly 
correlated to spin up (down) phases of the out of plane 
spin polarization. Since the autocorrelation function is 
related to the Fourier transform of the local density of 
states^, and because spin probes can be more demand- 
ing in practical implementations than charge detection, 
its evaluation yields useful indirect information on the 
spin degree of polarization. We believe these findings 
have the potential to provide interesting new strategies 
to dynamically control spin properties of charge carriers 
in graphene for future spintronics applications. 

The paper is organized as follows. In section II we de- 
scribe the spectum for static RSOC and introduce the 
model hamiltonian for periodically driven RSOC. Here 
we also present the exact solution to the dynamical equa- 
tions corresponding to the Dirac point k = 0. The main 
results of the Floquet-Magnus approach for the semi- 
analytical solution of the evolution operator at finite mo- 
mentum are presented in section III. Next, in section IV 
we compare the quasi-energy spectrum obtained through 
Magnus approach to that given by making a rotating 
wave approximation. We also evaluate and discuss the 
out of plane spin polarization as well as the onset of 
quantum revivals for the autocorrelation function. Fi- 
nally, in section V we give some concluding remarks and 
discuss an experimental scenario where our results could 
be tested. 



II. MODEL 

We consider a graphene monolayer sample subject 
to periodic time dependent spin-orbit interation of the 
Rashba type. In graphene RSOC interaction emerges as 
a consequence of a and 7r orbital mixing^ and stems 
from the induced electric field due to the substrate over 
which the graphene sample lies or by applied gate volt- 



ages. Then, a periodic modulation can in principle be 
implemented by means of time dependent gate voltages 
or by the induced time varying electric field within a par- 
allel plates capacitor coupled to a LC circuit. Under these 
circumstances the induced RSOC perturbation could be 
given a periodic time dependence V(t) = X(t)s. (z x a), 
where the driving amplitude will be assumed to be pe- 
riodic X(t + T) = X(t) with A(0) = \p, the coupling 
strenght in the static case. 

Concerning energy scales the value of the intrinsic and 
extrinsic spin-orbit coupling parameters A a nd X r in 
graphene have been obtai ned by tight binding^SI anc j 
band structure calculations^ 1 !. They gave estimates in 
the range 10 -6 — 10 _5 eV, much smaller than any other 
energy scale in the problem (kinetic, interaction and dis- 
order). However, the RSOC strength has recently been 
reported^ to be of order Xr rs 0.2t, where r rj 2.8eV 
is the value of the first-neighbor hopping parameter for 
graphene within a tight binding approach. 

The formulation of the problem is as follows. In mo- 
mentum space and and taking into account the energy 
scales of spin-orbit coupling we can work within the 
so called long wavelength approximation, where the to- 
tal hamiltonian for monolayer graphene in presence of 
time dependent RSOC can be described by the 8x8 
hamiltoniarP 

H(k,t) = (a x r z k x + <JyT y )s + X{t)(a x T z s y - a y s x ), (1) 

where vp ~ 10 6 m/s is the Fermi velocity in graphene, 
— {&x,<Jy,&z) is a vector of Pauli matrices, with 
(j z = ±1 describing states on the sublattice A(B) and so 
called pseudo spin degree of freedom, whereas r z = ±1 
describes the so called Dirac points K and K', respec- 
tively. In addition, k — (k x ,k y ) is the momentum mea- 
sured from the K point and Si (i=0,x,y,z) represents the 
real spin degree of freedom, with sq the identity matrix. 
In addition, A(t) gives the time dependence of the RSOC 
and we have neglected the intrinsic spin-orbit contribu- 
tions. Now since RSOC does not mix the valleys, we can 
focus on any of the two Dirac points, say K and then 
the results for the K' point are found by the substitu- 
tion k x — > —k x . Yet, we will formulate the problem in 
an isotropic way such that the results for the K' Dirac 
point will inmediately follow. Before dealing with the 
time dependent problem we summarize the main results 
for static RSOC. 

The spectrum of the noninteracting Hamiltonian 

Hq = hvpa . k (2) 
is given by the linear dispersion relation 



e°(fc) = ahvp^Jk^. + k^ = ahvpk, (3) 
whereas its eigenbasis is spanned by the spinors 
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with tan# = k y /k x , and a = 1 (—1) describes electron 
(hole) states. 

When a static RSOC interaction term is present the 
hamiltonian near the K point reads 

H(k) = hv F (a x k x + a y ky)s Q + Xr{cj x s v - a y s x ) (5) 
and the energy spectrum changes to ±e± with 

(6) 



e±(fc) = ±X R + JX 2 R + (Hv F k) 



Since RSOC mixes the a and 7r atomic orbitals it in- 
duces a gap Sq = 2A at the Dirac point k = 0. The 
static Rashba hamiltonian in Eq.([5| is diagonalized by 
the unitary transformation U(k) given explicitly as 



mi) = ± 



— ie i9 sin7 + — cos 7+ zcos7+ e~ j6, sin7 + \ 

— ie ie sin7_ cos7_ — icos7_ e -l9 sin7_ 

ie l6l sin7_ — cos7_ — icos7_ e~ l0 sin7. 
y ie ie sin7+ cos7 + - ~~' r : " ~ 



(7) 



i cos 7+ e %B sin 7+ J 



where cos j± = e± / J (hv F k) 2 + 1 
static RSOC hamiltonian reads 



In this basis the 



H(k) = Diag{-e + ,e_,-e_,e + }. 



(8) 



This particular choice of basis will simplify the calcula- 
tions that follow. 

We are interested in analyzing the emergent dynam- 
ics of Dirac fermions in monolayer graphene when the 
amplitude of RSOC is a periodically varying function of 
time X(t) = X R cos(Qt), with A^ and f2 the amplitude 
and frequency of the driving term. Then we would have 
to deal with the following 4x4 evolution equations 

ihdtVQc, t) = H(k, t) (9) 
However, if we make use of the unitary transformation 



([7]) the time dependent hamiltonian (JlJ becomes isotropic 
and block-diagonal 



H(k,t) 



h_(k,t) 






-(M) 



(10) 



with both sub-blocks periodic functions of time i.e. 
h±(k,t + T) = h±(k,t). Therefore the unitary trans- 
formation represented by ^ simplifies considerably the 
mathematical resolution of the evolution equations by re- 
casting the problem as two time reversal pairs of coupled 
2x2 two level problems. In addition, it has the physical 
appealing feature of clearly giving the subset of states 
which are dynamically coupled through the time depen- 
dent interaction. 

Let us then focus on the upper block h- (k, t) which 
reads 



h-(k,t) 



(hv F k) 2 + A(*)e+ hv F k[X R - X(t)} 
hv F k[X R - X{t)} ~(hv F k) 2 + X(t)e- 



(11) 



whereas the lower block is obtained by changing the sign 
of the amplitude X R — >■ —X R . Because of this symme- 
try relation among the two subspaces their quasi-energy 
spectra are identical. This is to be expected since RSOC 
is time reversal invariant. 

First of all, we notice that in the static limit A(t) — > X R 
the reduced hamiltonian ( 10 1 is block diagonal 



M*0 



(12) 



We also note that at the Dirac point k = 0, one gets 



MM 



-2A(t) 




In this case the resulting dynamics 

»ftfl t |0(f)> = Mo,t)Wt)) 

is exactly solved by the eigenspinors 
\Mt)) = (e M '« 0) 

\M*)) = (0,1) 



(13) 

(14) 
(15) 
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where 



/(*) = 



dt'A(t'). 



(16) 



As will be discussed below, these solutions correspond 
to zero quasi-energy Floquet states. The corresponding 
evolution operator is diagonal and given as 

[/_(0, t) = e l/( *)Diag{e l/(t) , e- l/(t) }. (17) 



III. MAGNUS-FLOQUET APPROACH 

Although we have shown that at k = the dynam- 
ics is exactly solvable, this is no longer true for finite k. 
Then, we need to resort to approximate solutions. As 
we discuss below, a semi-analytical approach known as 
Magnus-Floquet expansion will be suitable for dealing 
with the dynamical equations of periodically driven sys- 
tems. Since the Magnus-Floquet approach is not so pop- 
ular in the literature we now briefly summarize its main 
results (see reference^ for more detailed derivations). 

In the language of differential equations, the matrix so- 
lution S{t) of a n-dimensional system of dynamical evo- 
lution equations (here we omit the orbital degrees of free- 



dom for ease of notation), 

dMt) = A(t)9(t) 
i.e. a matrix that satisfies 

d t S(t) = A(t)S(t) 



(18) 



(19) 



is called a fundamental matrix solution if all its columns 
are linearly independent. If in addition, there is a time 
t = to such that 5 (to) is the identity matrix, then S(t) is 
call ed a principal fundamental matrix solution. To solve 
Eq.(19) Magnus 17 proposed to find exponential solutions 



to t 



le evolution operator in the form 



S(t) 



and then wrote M(t) as an infinite series 



M(t) = J2M j (t), 



(20) 



(21) 



where each term Mj (f ) is given as a combination of nested 
commutators, with the first terms reading as 



Afi(t) 
M 2 (t) 
M 3 (i) 



A(ti)dt x 



dh / [A(ti),A(*a)]*a 



6 Jo 



dti 



dt 2 / ([^o.^teMfo)] 

Jo 



(22) 
(23) 
(24) 



On the other hand, for periodic driving^LLiJ A(t + T) — 
A(t), Floquet's theorem states that the principal funda- 
mental solution of the dynamical equations can be writ- 
ten as 

S(t) = P{t)e tF (25) 

where P and F are n x n matrices, such that P(t) is pe- 
riodic P(t + T) = P{t) and F is time independent. Flo- 
quet's theorem is the time dependent analog of Bloch's 
theorem in solid state physics for spatially periodic struc- 
tures and it provides a time dependent transformation 
such that the so called Floquet states evolve according 
to the time independent matrix F. This time dependent 
transformation is implemented by P(t). 

One important remark is in order since although the 
interaction A(t) is periodic, the corresponding evolution 
matrix S(t) is not, i.e. S(t + T) ^ S{t). In fact S(T) 



carries indeed non trivial information on the dynamics of 
the periodic system. The eigenvalues of F are called Flo- 
quet exponents p. These Floquet exponents can be found 
by diagonalizing S(T) = e TF . Yet, they are not uniquely 
defined since p — > p + Him IT leaves S(T) invariant. 

In order to determine those exponents one standard 
approach consists of performing an expansion in the (in- 
finite) eigenbasis of time periodic functions £jv(i) = e lNnt 
(Fourier modes). In that periodic basis the evolution op- 
erator is diagonalized and the Floquet exponents q n are 
the logarithms of the eigenvalues of the evolution oper- 
ator evaluated at t = T, i.e. S(T). Then, in order to 
deal with the infinite eigenvalue problem one resorts to a 
truncation procedure in order to determine the Floquet 
exponents. 

The Magnus approach avoids the need to solve the infi- 
nite dimensional eigenvalue problem and has the physical 
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virtue of preserving unitarity of the evolved state to any 
order in the expansion. The connection between Magnus 
expansion and Floquet's theorem is found in reference^ 
where the authors present a solution of the evolution 
equations that consists of writing the periodic part P(t) 
as an exponential 



P(t) 



= e n(t) 



n(t + T) = n{t), 



(26) 



and then they proceed to expand both operators Q (t) 
and F in power series 



fi(i) 



OO 

E 



F 



E^- 

3=1 



(27) 



Now, since S(t) is by construction a principal fundamen- 
tal matrix solution P(T) = P(0) is the identity matrix 
and one gets for all values of j 



F,=M 3 (T)/T. 



(28) 



Introducing the Bernoulli numbers Bi (Bq = 1, B\ = 
-1/2, B 2 = 1/6, Bi = -1/30...) such that B 2m +i = 
for to > 1, the exponent operator term contributions 
Qj (t) satisfy a recurrence relation in terms of two auxil- 
iary time dependent operators W(t) and T(t), according 
to the relations 



W*) = E y ( W P (*) + (-i) m ?f (t)) 



(n > 1), 
(29) 



In turn, the W's and T's are given through the iterative 
relations 



j 


3-1 

= E 

m=l 




(i < / < i - 1) 


(30) 


mil) 
j 


3-1 

= E 

m=l 


ro T^- 1 )! 

'%i j-m J 


(i < i < j - 1) 


(31) 




= A, 




(j > i) 


(32) 


T (o) 
j 


= ^, 


(j > o). 




(33) 



In practical calculations, the relations in the last two lines 
serve to initialize the iterative procedure. 



IV. RESULTS AND DISCUSSION 

In order to apply the previous results to our prob- 
lem one just has to make the following identifications 
S(t) -> U-(k,t), A(t) -> -ih-(k,t)/h. Then, the 
characteristic exponents are proportional to the quasi- 
energies p = —iq n . For our calculation we choose a time 
dependence of the form X(t) = \r cos fit, where f2 is the 
frequency of the driving and Ai? the RSOC strength and 
proceeded to evaluate Fj by the iteration procedure de- 
scribed in the previous section. 




FIG. 1: (Color online) First Brillouin zone for the quasi- 
energy spectrum as function of adimensional non interacting 
quasi particle energy. Due to the dynamical interlevel mixing 
the static gap So gets closed (colored, thick continous lines) 
as compared to the static interacting spectrum (gray, thin, 
dashed lines). Colored arrows depict the limit A — > where 
the highly oscillatory contributions tend to cancel and the non 
interacting spectrum is recovered (black, thin, dashed lines). 
The Floquet Fourier solutions (colored, thick, dashed lines) 
show qualitative agreement with the Magnus result however 
they converge slower for small A. We have expressed all quan- 
tities in units of r = 2.8eV, the first neighbor hopping param- 
eter within a tight-binding approach and set fi = 1. 



We have found a characteristic behavior of the (adi- 
mensional) quasi-energies e± = q±/Ml, which qualita- 
tively do not change when one goes beyond third order 
in the Magnus-Floquet expansion (in the appendix we 
briefly describe the calculations up to fifth order). For 
finite n they explicitly read as 

e± = ±v/k 2 (16k 2 A 2 + A 4 - 2A 2 + 1), (34) 

where we have defined the adimensional quantities k — 
vpk/n and A = \r/MI. 

In FIG.l we show the typical behavior of the quasi- 
energies e± (blue and red, thick continous lines) given 
in equations ( |34[ ) as function of the adimensional quasi- 
particles non interacting energy e°. We also show (gray, 
thin, dashed lines) the corresponding static eigenvalues 
e± of the RSOC hamiltonian as described by h-(k), as 
well as the eigenenergies of the non interacting hamilto- 
nian (thin, dashed, black lines). We have also included 
(red and blue, thick, dashed lines) the result from a 20- 
mode Fourier expansion of the quasi-energies. We have 
expressed all quantities in units of the first neighbor hop- 
ping parameter t = 2.8eV and for finite ft have set ft = 1 
and change the effective coupling A. This is true for all 
the figures within the dynamical case. 

Since h-(k,t) mixes the static eigenstates of h- (k) the 
relative phases among them are dynamically modulated 
giving rise to interference phenomena and we find that 
this leads to a dynamical closing of the original gap Sq. 



G 



Therefore, the exact solutions at k = correspond to 
vanishing quasi-energies (modulo hfl). We further find 
(see arrows in FIG.l) that in the limit A — > 0, correspond- 
ing to a highly oscillating field, e± — > e± because then the 
influence of the driving quickly tends to vanish on average 
and intuitively one expects to recover the non interact- 
ing linear spectrum. The result from a 20— mode Fourier 
expansion shows qualitative agreement to the Magnus- 
Floquct approach, however a small discrepancy is found 
and the Magnus result converges faster in the limit of 
highly oscillating fields A — > 0, as depicted by thick ar- 



rows in FIG.l. 

Within the Magnus-Floquet approach the evolution 
operator U-(n,t) corresponding to the hamiltonian 
(k, t) is found to be given as 

U- (k, t) = e i /( t ) e ^(«.*)^ e «(«.*)-P' (35) 

where p is a vector of Pauli matrices and the non vanish- 
ing components of the vector V are given as 



V x {K,i) 

Vy(K,t) 

V z (k, t) 



-kA sin fit 



y/n 2 + A 2 
2kA(1 - cos fit) 
A 2 sinl7i , 2 

OK 



(4k 2 - 2A 2 + 1 + A 2 cos fit) 



\/k 2 + A 2 



k 2 cos 



fit) 



(36) 
(37) 
(38) 



whereas those of v are 

v x {n,t) = 
v z (n,t) = 



tnA 



Vk 2 + A 2 

tn 2 
\/k 2 + A 2 



(l + 4k 2 - A 2 ) 
(1 - 5A 2 ) . 



(39) 
(40) 



At the Dirac point n — we find that v x ,v z ,V x ,V y 



all 



vanish, whereas V z = A sin fit = f(t), and thus we re- 



cover the exact solution (17). 



Before we proceed to evaluate other physical quanti- 
ties of interest we make a brief disgression on the impor- 
tance of taking into account the full time dependence 
of the RSOC. In particular, due to smallness of the 
RSOC strength one can try a rotating wave approxima- 
tion (RWA) where for a given finite value of k only near to 
resonance (fl ~ 2vpk) terms are kept in the interacting 
hamiltonian. Now we show that for the present model 
this approach gives unphysical results, for instance, it 
predicts a gap opening at k = 0. This in turn would im- 
ply finite quasi-energy mode contributions at the Dirac 
point which contradicts the previously described exact 
result. 

Since the results are easily found within the 4- 
dimensional formulation we briefly return to the original 
basis and make a full dimensional discussion. Within a 
RWA approach, the hamiltonian reads (here we omit the 
spatial degrees of freedom for ease of notation 



H 



RWA 



(t) = Hq + iA(er+s_e 



(41) 



which for a given value of k describes near resonance 
fl ~ 2vpk spin and pseudo spin flipping processes and 
neglects the so called secular or counter-rotating terms 
that oscillate rapidly. In this case the solution is exact 



and the adimcnsional quasi-energy spectrum reads 



.RWA 



= 2 V S res + 5<r(«0 



(42) 



where S res = 2k — 1 describes the resonance and g a (k) = 
A 2 -k + 2<j v / k(A 2 + 1) + A 4 . When we evaluate at k = 
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FIG. 2: (Color online) Quasi-energies as function of momenta. 
Gray continuos lines correspond to the static limit. The solu- 
tion within this approach describes an unphysical gap opened 
at k = (see main text). 



and finite fl one would get the gaps A = v/(A 2 + 1/4 for 
the originally gapped states and a dynamically opened 
gap on the static degenerate states A dyn = 1. These re- 
sults are in disagreement with the exact solution to the 
full equations ( 15 ) where we had found zero energy solu- 



tions at k = so this approach is not suitable to describe 
the dynamical features of the periodic driving. In FIG. 2 



we depict the quasi-energy spectrum as a function of non- 
interacting energies e° for this RWA solution. 

After this brief discussion we turn back our attention 
to the Magnus-Floquet solution in order to get addi- 
tional information on the dynamical behavior induced 
on the system. This can be found by evaluating some 
other physical quantities of interest from an experimen- 
tal point of view. First we analyse the out of plane 
spin polarization S z (k,t) = {^{k,t)\S z \^(k,t)), where 
S z = h/2ao <g) s z and <7o the identity matrix in pseudo 
spin space. Using again the transformation ^ the out 
of plane spin polarization reads S z (k) = U(k)S z W (k) 
and is found to be isotropic and block anti-diagonal 



S z (k) = 

with s± given explicitly as 

h 











2Vk 2 + A 2 



(43) 



(44) 



with pi a vector of Pauli matrices and po the two di- 
mensional identity matrix. The anti-diagonal structure 
of S z in this basis reflects the fact that spin polarization 
is not conserved in presence of RSOC. Then we have to 




FIG. 3: (Color online) Density plot showing the out of plane 
spin polarization plotted against normalized time t = h/r for 
the exact solution at « = 0. The periodicity is inherited from 
the driving field. The vertical axis indicates the normalized 
adimensional field strength A and one gets up (+) down (— ) 
spin components phases equally separated by the zeroes of 

/(*)- 



evaluate 

S z (k,t) = (V(k,0)\uHk,t)S z (k)U(k,t)\*(k,0)), (45) 




2n 



4n 



FIG. 4: (Color online) Density plot showing the time behav- 
ior of the out of plane spin polarization for the semi- ana- 
lytical Magnus-Floquet solution at finite k = 1. Alternating 
maxima and minima appear due to interlevel mixing among 
different static eigenstates. This is also manifested by the lost 
of symmetry with respect to t = n and stems from the driven 
modulated phases since dynamical interlevel mixing leads to 
interference effects. 



for any initially prepared state |^/(/c,0)). Next we sep- 
arate explicitly the four spinor 0)) in upper and 
lower components as 



|*(fc,0)) 



1 

71 



(46) 



where ip± are normalized two dimensional spinors. After 
some algebra one gets for the spin polarization in terms 
of adimensional parameters 



S z (n,t) = K^le" 



-2i/(t) 



-iv-p 
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For a finite value of k we now choose the initial spinor 
configuration ip± = (±i, l)/y/2, in such way that the out 
of plane polarization vanishes (\ i S(k,t)\S z (k)\ i If(k,t)) = 
for a static RSOC. 

In figure FIG. 3 we show a density plot of the result- 
ing out of plane spin polarization for the exact solution 
k = 0. In this case, the only relevant parameter is the 
adimensional amplitude of the driving field A. As ex- 
pected, for A = the system remains unpolarized for all 
values of the adimensional time t (given in units of t/H). 
For finite values of the effective coupling an alternating 
pattern of spin phases (denoted as + and — represent- 
ing up and down, respectively) are seen to apper as time 
evolves. They are symmetrically distributed among the 
values t = mr where f(t) and thus the relative phases 
among the static RSOC eigenstates vanish. 
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FIG. 5: (Color online) Density plot showing the behavior 
of the absolute value of the autocorrelation function plotted 
against normalized time t = h/r and strength A for the exact 
solution at k — 0. As expected, only for large values of the 
interaction strength the evolved state departs considerably 
from the initial state configuration (see main text). 



However, as shown in figure FIG. 4, once K is finite this 
panorama qualitatively changes. In this case, the addi- 
tional interference due to level mixing induces a pattern 
of alternating maxima (+) and minima (— ) for t n = nn, 
n =G N. The reason for such a behavior is that for a 
given t n the evolution operator is given by e TF and thus 
the polarization maxima and minima S(n, T) = ±1/2 de- 
pend essentially on the quasienergy spectrum properties. 
Then, increasing A makes these alternating maxima and 
minima to get closer. When we move to the next period 
(corresponding to t — 4n in figure FIG. 4) the number of 
alternating maxima and minima is doubled, and so on. 
Therefore, dynamical coupling produces a non vanishing 
out of plane spin polarization with an oscillating time 
pattern resulting from the mixing of the static eigenstates 
such that the changing relative phases, modulated by the 
time dependent interaction prevent total destructive in- 
terference to happen. Therefore ac driven RSOC pro- 
vides a suitable means to dynamically control the degree 
of spin polarization and could in principle serve to gen- 
erate non vanishing and non trivial spin polarized phases 
in otherwise unpolarized states in monolayer graphcnc. 

In order to complement the just described physical pic- 
ture of the spin polarization scenario we evaluated the 
autocorrelation function A(H,t). This is given by the 
projection of the evolved state \^(k, t)} along a given (in 
principle arbitrary) initial spinor configuration \^(k, 0)), 



FIG. 6: (Color online) Density plot showing the behavior of 
the absolute value of the autocorrelation function against nor- 
malized time t — h/r and RSOC strength A for the Magnus- 
Floquet solution at K = 1. The phase interference effects 
previously discussed. Now the onset of recurrences at small 
A is a consequence of interlevel mixing. 



i.e. 

A(K,t) = <*(*> 0)|¥(/5,t)). (48) 

The absolute value of the autocorrelation function pro- 
vides information on the so called recurrences or quantum 
revivals of the dynamics, i.e. those values of the time pa- 
rameter for which such overlapping is a maximum. In 
addition, its Fourier transform is proportional to the lo- 
cal density of states^. 

In figure FIG. 5 (FIG. 6) we plot the absolute value of 
A(n, t) obtained by means of the exact (semi-analytical) 
evolution operators. For the exact solution shown in fig- 
ure FIG. 5 only large values of A induce an appreciable 
phase change and the system remains mostly correlated 
to the initial state. As for the case of the out of plane spin 
polarization, for t — t n the autocorrelation gives max- 
ima corresponding to red (dark gray) zones in the figure. 
These signal the return of the system to initial vanishing 
spin polarization. Maxima and minima of spin polar- 
ization correspond to partial quantum revivals. Given 
the definition of the autocorrelation as the probabibility 
that the system returns to its initial state we find that for 
A = giving e — k the loci of A(k, t n ) — 1 correspond to 
vanishing of spin polarization and are thus equally spaced 
as St = 7T. On the contrary, the maxima and minima of 
spin polarization correspond in this case to quantum sup- 
pressions and are shown as blue (dark) zones in the figure. 
As soon as we move towards finite values of k (FIG. 6) in- 
terference phenomena come again into play and we find 



recurrences represented as read (dark gray) zones and 
suppressions, described by purple (black) zones. These 
arise because of the constructive and destructive inter- 
ference effects, described previously and are modulated 
by the ac driving. To check that this physics is different 
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FIG. 7: (Color online) Density plot showing the behavior of 
the absolute value of the autocorrelation function versus nor- 
malized time and RSOC strength A for the static interaction 
at K = 0. In this case, quantum revivals are separated by 
the inverse of the gap So = 2A in accordance to Heisenberg 
uncertainty principle. 



from the static scenario, in Figures FIG. 7 and FIG. 8 we 
depict the corresponding countour plots for time inde- 
pendent RSOC interaction. For the k = case shown in 
FIG. 7, the alternating pattern of quantum revivals and 
suppressions is separated by the inverse of the energy gap 
Sq. This is due to the fact that the in this situation the 
relative phase among the non interacting eigenstates is 
set by the energy separation among them, i.e. So, and ac- 
cording to Heisenberg's time-energy uncertainty relation 
one should have St inversely proportional to the energy 
separation. As soon as k is finite, a quasi homogeneous 
pattern of recurrences is seen. Therefore although in the 
static case where the spin polarization vanishes for all 
values of RSOC strength, there appear recurrences indi- 
cating there is no longer correlation between S z and A(t), 
which for the system under study are only correlated for 
time dependent driving. 

Two comments are in order here. The first one is con- 
cerning the validity of the Dirac fermion hamiltonian ap- 
proximation in order to deal with spin-orbit related phe- 
nomena in graphene monolayer. As it is discussed in 
reference^ the large gap of the hamiltonian describing 
the a orbitals implies that the effective hamiltonian in- 
cluding SOC is essentially given by the long wavelength 
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FIG. 8: (Color online) For finite values of momentum k = 1 
recurrences are modulated signaling the coexistence of re- 
vivals for times inversely proportional to the static gap So = 
2A and those corresponding to the larger energy separation 
5 = \/k 2 + A 2 , which correspond to a smaller time and seen 
in the figure as modulated alternating pattern of maxima and 
minima for each recurrence curve corresponding to k = 0. 



or Dirac Hamiltonian considered in our model. In addi- 
tion, M. Gmitra, et aP^used first principles calculations 
to discuss the relevance of spin-orbit related physics in 
graphene near and beyond the Dirac points. The second 
point to remark is the role of localized impurities. As 
has been discussed irP3 the presence of local impurities 
can enhance the value of the static intrinsic spin-orbit 
coupling strength because of the induced sp 3 distortion 
which leads to a hybridization of the tt and a orbitals and 
as shown in reference^ RSOC only respond to this a — ir 
hybridization. Therefore, we would expect our results to 
be robust or even enhanced if localized impurities were 
included in the model. 

Now we would like to compare to other proposals of 
dynamical modulation of energy gaps under ac drivenin 
graphene. For instance, Oka and AokP^l found that cir- 
cularly polarized intense laser fields can induce a photo- 
voltaic effect in graphene, that is, a Hall effect without 
magnetic fields. This in turn relies on the gap opening 
at the Dirac point k = 0. However, as shown recently 
by Zhou and WtP^ in analyzing the optical response 
of graphene under intense THz fields, the contributions 
from large momenta make the effective gaps opened to 
get also dynamically closed. However, it is found that the 
quasi energy spectrum for the linear polarization leads 
to linear quasi energy spectrum (see alsrJ^). Our model 
could in principle be mapped into their scenario of lin- 
early polarized radiation field. Yet, in both workiP^ZI 
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the case of linear polarization leads to linear quasienergy 
spectrum. Therefore, from the semi parabolic quasi en- 
ergy spectrum shown in figure (FIG IV ) we can infer 
that the bending of the quasienergy spectrum makes 
the spin-orbit driven scenario qualitatively different from 
these other approaches to dynamical control of graphene 
electronic properties. Our argument relies on the fact 
that the topological properties of periodically driven sys- 
tems are characterized by the quasienergy spectrum^ so 
we can conclude that the physics related to ac driven 
spin-orbit does reveal new physical interesting electron- 
ics properties which are absent in the static regime. 



V. CONCLUSIONS 

We have investigated the role of periodically driven 
RSOC in monolayer graphene and recasted the original 
4-dimensional problem as an equivalent set of two two- 
level problems. Due to the induced modulation of the 
relative phases among the static eigenstates we found a 
closing of the static gap at the Dirac point k = 0. This 
result is in agreement with the available exact solution 
and differs from RWA where an unphysical gap is seen 
to appear. This physical picture is confirmed through a 
Fourier mode expansion and we found that Magnus Flo- 
quet approach indeed has the advantage of providing the 
quasienergy spectrum with less computational effort. We 
also found that the generation and manipulation of out 
of plane spin polarization for otherwise spin unpolarized 
states requires the time driving to be realizable within 
this set up. Due to the induced interlevel mixing among 
the static eigenstates we found a set of alternating posi- 
tive and negative spin phases in clear distinction to the 
well separated spin phases at the Dirac point correspond- 
ing to the exact solution. The dynamical onset of quan- 



tum revivals described through the autocorrelation func- 
tion is directly correlated to the appearance of maxima 
for either spin phases. However, in the static case, such 
a correlation does not ensues since the spin polarization 
vanishes identically whereas the quantum revivals are 
still present. Concerning the actual experimental realiza- 
tion of our proposal we believe it could be implemented 
by means of Magnetic Resonance Force Microscopy as 
reported irPH Within this scheme, single spin polariza- 
tion could be detected by means of the frequency shift 
induced on a cantilever that is used to scan the sample. 
The sign of the cantilever's frequency shift can be associ- 
ated to the spin polarization. This detection is achieved 
by means of low intensity magnetic fields under resonant 
conditions, thus no magnetic coupling terms need to be 
included in the description of the dynamics of the charge 
carriers in graphene. In this way, although in principle 
RSOC in graphene has a small static strength, the time 
dependent efective phenomena produce some interesting 
spin controlling strategies which we believe could provide 
a route to new implementations of graphene in spintronic 
devices with appropiate spin detection techniques. 
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Appendix A: derivation of iterative terms and quasi 
energy spectrum 

Using the simplifying notation for the time dependent 
interaction h~ (fc, tj) = one gets to fourth order in the 
iteration procedure that 



Fi 
F 2 
F 3 
Fi 



Th Jo 



—i 

2TH 2 
—i 



dti / /i^]di2 
o Jo 

3 rT 



{)Th , j dh J dt 2 J , [h 2 ,h 3 }} + [h 3 , [h 2 , h x ]])dt a 



dh / {[[[hx ,h 2 ],halK] + 
Jo 



dti dti 
12Th*J Vo Jo 



(Al) 

(A2) 
(A3) 
(A4) 

(A5) 



and performing the corresponding calculations one gets whereas the even contributions F 2 j all vanish. The 
. ^ quasienergies e± are obtained through diagonalization of 

F ^ - -jwrm {K(Tz + A(T * } (A6) 
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the selfadjoint matrix — iF = —i J2j Fj an d up to this 
order one gets 

e± = ±-*-[k 2 (331776k 6 A 2 + 81(A 2 -2) 4 - (A9) 
3d 

2304k 4 (23A 4 - 72 A 2 ) + 32k 2 (1109 A 6 
-900A 4 - 324A 2 )] 1/2 
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